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ABSTRACT 

Aims. We report on observations of the eclipsing and interacting binary /? Lyrae from the Suzaku X-ray telescope. This system involves 
an early B star embedded in an optically and geometrically thick disk that is siphoning atmospheric gases from a less massive late 
B II companion. 

Methods. Motivated by an unpublished X-ray spectrum from the Einstein X-ray telescope suggesting unusually hard emission, we 
obtained time with Suzaku for pointings at three different phases within a single orbit. 

Results. From the XIS detectors, the softer X-ray emission appears typical of an early-type star. What is surprising is the remarkably 
unchanging character of this emission, both in luminosity and in spectral shape, despite the highly asymmetric geometry of the system. 
We see no eclipse effect below 10 keV. The constancy of the soft emission is plausibly related to the wind of the embedded B star and 
Thomson scattering of X-rays in the system, although it might be due to extended shock structures arising near the accretion disk as 
a result of the unusually high mass-transfer rate. There is some evidence from the PIN instrument for hard emission in the 10-60 keV 
range. Follow-up observations with the RXTE satellite will confirm this preliminary detection. 
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1. Introduction 

The binary ft Lyr is a nearly edge-on, semi-detached interact- 
ing system that has undergone mass reversal and remains in a 
phase of large-scale mass transfer. The primary, mass-losing star 
(the "Loser") is a B6-B8 Hp star. The mass-gaining star (the 
"Gainer") is embedded in an optically thick accretion disk and 
is not directly visible. Although the embedded source had been 
considered as a possible compact object (Devinney 1 97 1 ; Wilson 
1971), it is probably a main sequence B0 star (Hubeny & Plavec 
1991). The system is very complex, having bipolar jet-like struc- 
tures (Harmanec etal. 1996; Hoffman et al. 1998), a circumbi- 
nary envelope (Batten & Sahade 1973; Hack et al. 1975), and a 
substantial kilo-Gauss magnetic field (Leone et al. 2003). 

The optical light curve of the system features a primary min- 
imum that is ss 1 magnitude deep and a secondary minimum 
« 0.4 magnitudes deep (see Fig. [TJ; however, the secondary 
minimum is deeper than the primary minimum at shorter wave- 
lengths, and below Lycr the eclipses no longer appear (Kondo 
et al. 1994). A summary of the system properties is given in 
Table Q] The orbital period is 12.9 days, and the mean light 
curves appear stable with epoch. The UV spectrum of ft Lyr is 
dominated by an anomalous continuum and emission lines with 
unusually strong P Cygni profiles typical of hot star winds (Hack 
et al. 1975; Aydin et al. 1988; Mazzali 1987). 
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Despite numerous and ongoing modeling attempts (e.g., 
Wilson 1974; Linnell & Hubeny 1996; Bisikalo et al. 2000; 
Linnell 2002; Nazarenko & Glazunova 2003, 2006ab), no model 
is yet capable of matching the observed light curves of ft Lyr 
from the IR through the UV. Strangely, ft Lyr has been largely 
unstudied in the X-ray regime, despite the strong X-ray flux de- 
tected by the ROSAT HRI (Berghofer & Schmitt 1994). An un- 
published spectrum taken with the EINSTEIN/SSS in 1979 re- 
veals X-ray emission at relatively high energies, suggesting that 
phase-dependent observations may provide new clues to resolv- 
ing the puzzle of the ft Lyr geometry and interactions. Exploiting 
Suzaku 's excellent sensitivity to hard X-rays, we conducted three 
pointed observations of ft Lyr within the same orbit. In the fol- 
lowing section the observations and reduction of data are de- 
tailed. Analyses of the spectra with phase are described in §3, 
and a discussion of the results is presented in §4. 



2. Observations and Data Reduction 

The joint Japan/US X-ray astronomy satellite Suzaku (Mitsuda 
et al. 120071 1 observed ft Lyr in May 2006 on three occasions at or- 
bital phases as shown in Fig.Q] with corresponding viewing per- 
spectives of the ftLyr system illustrated in Fig. [2] (Note that the 
primary minimum occurs when the loser star is eclipsed, and the 
secondary minimum when the disk component is eclipsed.) The 
pointings were spaced approximately 4.3 days apart to sample 
the full 12.9 day orbit of the binary. Exposure times and count 
rates are tabulated in Table [2] 
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Table 1. Properties of (3 Lyrae 



cp = 0.24 



Component" 


Gainer Loser 


r cff [K] 


32,000 13,300 


Spectral Type 


»B0 V B6-8 Hp 


M/M G 


^13 ~3 


M b [M a yr 1 ] 


4.7 x 10~ 8 7.2 x 10~ 7 


!^[kms-l] 


1470 390 


System Property" 


Value 


Orbital Period 


12.9 days 


Viewing Inclination 


SO 


Binary Separation 


55 - 60R o 


Distance 


270 pc 


\ogN H (cm -2 ) 


20.76 


ROSAT PSPO: 


0.07 cps 


Einstein SSS'': 


0.11 cps 



" Component and system properties taken from 
Harmanec (2002) unless otherwise noted. 
* Mass-loss rate and terminal speed for the stel- 
lar winds of the respective binary components 
(Mazzali 1987). 
c Berghofer & Schmitt 1994) 
d Waldron (private comm.) 
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Fig. 1. Illustration of the Suzaku pointings with respect to the 
binary phase. The solid line represents the normalized V-band 
light curve (Fourier fit by Harmanec et al. 1996). Primary eclipse 
(phases 0.0 and 1 .0) occurs when the mass-losing star is occulted 
by the disk. The three intervals marked on the curve represent 
our Suzaku pointings, which occurred during a single orbit in 
May 2006. Phases were calculated from the quadratic ephemeris 
of Harmanec & Scholz (1993). 



Suzaku carries four X-ray Imaging Spectrometers (XIS; 
Koyama et al. |2007| > and a collimated Hard X-ray Detector 
(HXD; Takahashi et al. |2007l >. The field-of-view (FOV) for the 
XIS detectors is 17' x 17'. One of the XIS detectors (XIS1) 
is back-side illuminated (BI) and the other three (XIS0, XIS2, 
and XIS3) are front-side illuminated (FI). The bandpasses are 
-0.4- 12keV for the FI detectors and -0.2- 12keV for the BI 
detector. The BI CCD has higher effective area at low energies, 
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Fig. 2. Graphic topview representation of the fi Lyr star-disk sys- 
tem, following the model of Hoffman et al. (1998). Components 
are illustrative and not to scale. The three Suzaku pointings oc- 
curred at phases ip = 0.24, 0.55, and 0.91 all within the same 
orbit. Arrows indicate the Earth's line-of-sight for these three 
phases. 



however its background level across the entire bandpass is higher 
compared to the FI CCDs. 

The angular resolution of the X-ray telescope onboard 
Suzaku is w 2'. Therefore in the XIS image, /?Lyr is not re- 
solved from two nearby B-type stars: HD 174664 (jSLyrB) and 
HD 174639. While the latter star was not detected by ROSAT, 
the former has a ROSAT HRl count rate of 4 x 10~ 3 cps, as com- 
pared to the /?Lyr ROSAT HPJ count rate 4 x 10" 2 cps. We are 
confident that the X-ray flux detected by Suzaku is at least 90% 
dominated by fi Lyr. 

The HXD consists of two non-imaging instruments (the PIN 
and GSO; see Takahashi et al. 120071 1 with bandpasses of ~ 10- 
70 keV (PIN) and ~ 40-600 keV (GSO), and a FOV of 34' x 34' 
(PIN). Both of the HXD instruments are background-limited. 
The background subtraction for the HXD is performed by mod- 
eling the background spectrum. Presently, the non-X-ray back- 
ground model (e.g., particle events) is known for the PIN de- 
tector with ~ 3-5% accuracy (Kokubun et al. |2007| l. The back- 
ground modeling for GSO data is currently far less certain than 
for the PIN, and so we do not report on the GSO measurements. 

HXD-PIN data reduction and extraction of spectra were 
performed using the latest calibration sources and background 
models. We account for cosmic X-ray background (CXB) in 
fitting the spectral models using the procedure suggested by 
the Suzaku team based on a "typical" CXB spectrum (see 
www.astro.isas.ac.jp/suzaku). We do not correct for contribu- 
tions to the HXD background due to the South Atlantic Anomaly 
(SAA; see Kokubun et al. |2007t because our observations of 
jSLyr were performed when the background count-rate was at 
its lowest. 

Based on the Rosat All-Sky Survey (RASS), there are a few 
X-ray sources in the 34' x 34' PIN's field of view (FOV); how- 
ever, our target ySLyr is by far the brightest. The stellar coro- 
nal X-ray sources present in the FOV are expected to be faint 
in the HXD's energy range. There is also an active galactic nu- 
cleus (AGN) in the FOV about 15' away from fi Lyr. This AGN, 
QSO B 1847+3330, is cataloged in the RASS with a count-rate of 
0.05 cps. To estimate its potential contribution to the 10-70 keV 
energy range, we adopted a standard AGN power-law spec- 
trum with F = 2 and a low interstellar absorption column of 
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Date 




Exposure 




XIS Count rate (10 2 cps) 








(ksec) 


XIS0 


XIS1 


XIS2 


XIS 3 


2006 May 7 


0.55 


15.4 


6.35 ± 0.24 


8.53 ±0.36 


6.33 ± 0.24 


5.15 ±0.22 


2006 May 12 


0.91 


17.7 


5.76 ±0.21 


8.14 ±0.33 


5.56 ±0.21 


4.86 ± 0.20 


2006 May 16 


0.24 


15.7 


6.15 ±0.23 


9.58 ±0.36 


6.23 ± 0.22 


4.93 ± 0.22 



" The orbital phases for /? Lyr were computed at the midpoint of each observation using the quadratic 



ephemeris of Harmanec & Scholz (1993). Phases 0.0 and 1.0 correspond to primary eclipse. 



Table 3. PIN Observations of j3 Lyrae 





Exposure 


Count Rate 


Background 




(ksec) 


(cps) 


(cps) 


0.55 


14.7 


0.623 ± 0.007 


0.536 ±0.002 


0.91 


17.0 


0.505 ± 0.005 


0.509 ± 0.002 


0.24 


14.1 


0.573 ± 0.006 


0.564 ± 0.002 



Nn = 5 x 10 20 cm 1 . The predicted count-rate for the Suzaku PIN 
is « 0.006 cps. The observed HXD-PIN count-rates are listed in 
Tabled The contribution of QSOB 1847+3330 is between 10- 
30% of the detected flux. 



3. Analysis 

A comparison of our three XIS spectra below 10 keV with the 
XIS detectors indicates very little variability with phase (see 
Fig. [3}. Based on chi-square model fits to the data, there is no ev- 
idence for statistically significant differences between the three 
spectra. For illustration a model fit for the ip = 0.24 pointing 
appears as the solid line in all three panels of Fig. [3j its good 
agreement with the data at all phases demonstrates that the soft 
X-ray spectrum of /3 Lyr is nearly constant in shape and strength. 
Independent fits for each pointing are very similar. 

These XIS X-ray spectra are most probably thermal in na- 
ture, since emission lines are detected at 1.35 keV (Mg xi) and 
1.86 keV (Si xm). Our fits indicate a two-component model with 
the majority of the XIS X-ray emission arising from a tempera- 
ture of as 7.2 MK and a hotter but much weaker component of 
> 20 MK. Solar abundances were adopted except for nitrogen 
that was enhanced by more than 10 times to achieve an ade- 
quate fit. The derived hydrogen column density from the model 
is Nh — 6.5 + 0.24 x 10 20 cm -2 , largely consistent with the in- 
terstellar value. 

Hard emission in the 10-60 keV band of the PIN is back- 
ground dominated. The observed total and the background count 
rates are given in Table [3] The source counts come from the dif- 
ference of the total and background values. Based on the total 
counts, there is a significant source detection at phase <p = 0.55, a 
non-detection for tp = 0.9 1 , and a marginal detection at <p = 0.24. 
Emission at such hard energies is exceptionally unusual and 
unexpected for a system like /3 Lyr that consists of two early- 
type stars. The implied hard X-ray luminosity above 10 keV at 
<p ~ 0.55 is about 10~ 3 L Q . If confirmed, this hard component 
would be an important tracer of the flow geometry and shock 
structure since the circumstellar material is quite transparent to 
such high energy X-rays. 




Energy (keV) 



Fig. 3. Spectra of fl Lyr from the XIS 1 detector for each pointing. 
Phase within the orbit is indicated in each panel. The data are of 
reasonably high signal-to-noise, and it is clear that the spectral 
distribution shows little change with phase, being of nearly con- 
stant brightness and spectral shape. The solid line is a model fit 
to the data of <p — 0.24 that is replotted in the other phases for 
comparison. 

4. Discussion 

The results of our Suzaku study are both revealing and perplex- 
ing. We had expected to see an eclipse of soft X-rays and es- 
sentially no hard emission. Instead, we found that the soft X- 
rays were nearly constant, and there is an indication of a quite 
hard emission component in the system. The near constancy of 
the soft spectrum suggests that the X-ray source must be axially 
symmetric. This is quite surprising, since the system is intrinsi- 
cally non-axisymmetric. 

It is worth recalling that the UV lightcurve below Lya is 
notable for lacking any eclipses as well (Kondo et al. 1994). A 
somewhat analogous source to /3 Lyr is the near edge-on interact- 
ing binary W Ser. Weiland et al. (1995) interpret UV spectra of 
W Ser in terms of an extended boundary layer between the star 
and the accretion disk. Thus, it may be possible that the largely 
steady FUV and soft X-ray emission could arise from a similar 
region in the /3 Lyr system. However, we do not currently favor 
this explanation since the most recent modeling of the disk and 
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star components of fi Lyr (Linnell 2002) indicate that gainer star 
is entirely obscured at all phases. 

An alternative view is suggested by hydrodynamic simula- 
tions. Bisikalo et al. (2000) modeled the mass transfer of fi Lyr 
and found that they could reproduce a disk-like structure but 
without a well-identified "hot spot". Instead, their simulations 
yielded a complex distribution of extended shocks. They also 
were able to generate a bipolar flow similar to the jet detected 
by Harmanec et al. (1996) and Hoffman et al. (1998). However, 
radiative cooling was not treated self-consistently in those sim- 
ulations, but approximated through the equation of state. It is 
unclear whether this model can generate gas at sufficiently high 
temperatures to produce the observed X-rays. Hydrodynamic 
modeling of the fi Lyr system was also pursued by Nazarenko 
& Glazunova (2003; 2006ab), first in 2D simulations and then 
in 3D simulations. These researchers included cooling curves in 
their simulations and allowed for a stellar wind by the gainer. 
Their models also lead to a disk, a bipolar flow, and an environ- 
ment permeated with elongated shocks at different azimuthal ori- 
entations about the gainer. "Warm" gas temperatures up to about 
200,000 K are achieved, but such gas can not contribute signifi- 
cant emission to the X-ray band. The gas dynamical simulations 
are qualitatively promising in terms of predicting shocked struc- 
tures that are spatially distributed in radius and azimuth through- 
out the disk. In such a model, the hot gas could potentially that 
hot gas could be viewable at every phase. However, X-ray emis- 
sions were not the focus of those studies, and it would be useful 
to have new simulations that emphasize the hot plasma struc- 
tures. 

We prefer an interpretation in which the X-rays arise 
from distributed shocks in the wind of the gainer star. At 
kT « 0.6 keV, the spectral characteristics of the Suzaku spectra 
are compatible with a typical early main sequence B star, and the 
XIS luminosity Lx(XIS) « 6.6 x 10 30 erg s _1 is commensurate 
with the soft emission expected from an early B0-B1 star 
(Cohen, Cassinelli, & MacFarlane 1997). Although the favored 
model for the binary geometry may preclude a direct view of 
the gainer star, it does permit a reasonably deep view into its 
wind, even at secondary eclipse (with the giant star at front). 
Moreover, we know from the polarimetric study of Hoffman 
et al. (1998) that there is substantial scattering opacity above 
and below the disk plane. The near constancy of the X-rays seen 
with the XIS may result from an extended "halo" of scattered 
X-ray light, similar in spirit to the scattered optical light 
observed in images of some Herbig-Haro objects seen edge-on 
to their disks, such as HH 30 (Burrows et al. 1996). In this 
picture the soft X-rays ultimately originate in the shocked wind 
of the early gainer star. Part of the emission is observed directly 
from hot plasma at large radii in the wind, and part is scattered 
into our line-of-sight from above and below the disk. Thus the 
0.6 keV temperature of the hot gas would represent an upper 
limit value owing to the fact that photoabsorption of X-rays 
is more severe toward the soft end of the spectrum. On the 
other hand, the observed X-ray luminosity must be a lower limit. 

We have considered three distinct origins for the observed 
X-ray emissions. At this point we favor the gainer wind as the 
source of the soft X-rays, but acknowledge that future hydro- 
dynamic modeling may change this picture. In addition, there 
seems to be an unusually hard component of emission in the sys- 
tem, in excess of 10 keV. Such hard emission is more typically 
associated with compact objects, and not the winds of hot stars 
(e.g., Berghofer et al. 1997). However, the nature of the gainer 
is still ambiguous; perhaps the X-ray properties of /?Lyr could 



be understood in relation to a central compact object (Devinney 
197 1 ; Wilson 197 1). Although this seems unlikely, the X-ray ob- 
servations certainly do not conform to our original expectations, 
and so at least a reconsideration of the possibility seems in or- 
der. The X-ray pointings by Suzaku offer the tantalizing promise 
of providing new and valuable information about the fi Lyr sys- 
tem, but it is clear that a far more rigorous sampling of the X-ray 
lightcurve and source spectrum and new models for the system 
will be needed to test the suppositions that we have put forth. 
We have obtained time with the RXTE satellite to create a more 
complete X-ray light curve of fi Lyr, which should provide new 
insight into this complex system in the near future. 
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